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Thermoelectric  generator  (TEG)  has  attracted  considerable  attention  for  the  waste  heat  recovery  of  inter¬ 
nal  combustion  engine.  In  this  study,  a  3-D  numerical  model  for  engine  exhaust-based  thermoelectric 
generator  (ETEG)  system  is  developed.  By  considering  the  detailed  geometry  of  thermoelectric  generator 
(TEG)  and  exhaust  channel,  the  various  transport  phenomena  are  investigated,  and  design  optimization 
suggestions  are  given.  It  is  found  that  the  exhaust  channel  size  needs  to  be  moderate  to  balance  the  heat 
transfer  to  TEG  modules  and  pressure  drop  along  channel.  Increasing  the  number  of  exhaust  channels 
may  improve  the  performance,  however,  since  more  space  and  TEG  modules  are  needed,  the  system  size 
and  cost  need  to  be  considered  as  well.  Although  only  placing  bafflers  at  the  channel  inlet  could  increase 
the  heat  transfer  coefficient  for  the  whole  channel,  the  near  wall  temperature  downstream  might 
decrease  significantly,  leading  to  performance  degradation  of  the  TEG  modules  downstream.  To  ensure 
effective  utilization  of  hot  exhaust  gas,  the  baffler  angle  needs  to  be  sufficiently  large,  especially  for 
the  downstream  locations.  Since  larger  baffler  angles  increase  the  pressure  drop  significantly,  it  is  sug¬ 
gested  that  variable  baffler  angles,  with  the  angle  increasing  along  the  flow  direction,  might  be  a  middle 
course  for  balancing  the  heat  transfer  and  pressure  drop.  A  single  ETEG  design  may  not  be  suitable  to  all 
the  engine  operating  conditions,  and  making  the  number  of  exhaust  channels  and  baffler  angle  adjustable 
according  to  different  engine  operating  conditions  might  improve  the  performance. 

©  2014  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

In  recent  years,  the  global  energy  and  environmental  issues 
have  driven  the  application  of  waste  heat  recovery  techniques  in 
main  energy  conversion  devices  such  as  internal  combustion 
engine  (ICE).  For  a  typical  gasoline  engine  as  an  example,  only 
25%  of  the  chemical  energy  in  fuel  can  be  converted  into  useful 
energy  for  powering  vehicles  and  other  devices,  and  the  rest  of 
the  energy  is  discharged  to  the  ambient  environment  in  the  form 
of  waste  heat  through  the  exhaust  and  coolant  [  1  ].  If  6%  of  the 
exhaust  heat  is  converted  into  electrical  energy,  it  would  be 
possible  to  reduce  the  fuel  consumption  by  about  10%  [2],  Thus, 
the  technologies  of  waste  heat  recovery  for  ICE  become  signifi¬ 
cantly  important  for  improving  the  overall  efficiency  and  reducing 
the  fuel  consumption.  Among  the  different  technologies,  thermo¬ 
electric  generator  (TEG)  has  received  many  attentions  due  to  the 
advantages  such  as  free  maintenance,  high  reliability  and  quiet 
operation. 
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Many  studies  have  been  conducted  to  use  TEG  for  the  waste 
heat  recovery  of  ICE  exhaust.  LaGrandeur  et  al.  [3]  developed  and 
tested  a  TEG  system  based  on  an  in-line  six-cylinder  gasoline 
engine  in  BMW  5301.  It  was  found  that  the  influence  of  the  primary 
heat  exchanger  backpressure  should  not  be  neglected  in  high-load 
operation  (e.g.  highway  driving  cycles).  Thacher  et  al.  [4] 
conducted  tests  of  an  automobile  exhaust  TEG  system  installed 
in  a  CMC  Sierra  pickup.  They  concluded  that  the  performance  of 
the  system  is  sensitive  to  the  coolant  temperature,  and  the  system 
weight  should  be  considered.  Similarly,  Hsu  et  al.  [5]  tested  the 
performance  of  a  Bi2Te3-based  TEG  system,  which  is  connected  to 
the  middle  of  an  exhaust  pipe  in  a  Chrysler  Neon.  Saqr  et  al.  [6] 
made  a  comprehensive  review  of  the  main  aspects  of  exhaust- 
based  thermoelectric  generator  (ETEG)  in  the  past  20  years,  and 
summarized  the  main  challenges  for  ETEG  systems  such  as  the 
need  in  efficiency  improvement. 

Beside  the  above  mentioned  tests  of  ETEG  systems  in  vehicles, 
in-depth  investigation  and  optimization  of  ETEG  system  are 
necessary  for  performance  improvement.  Typically,  an  ETEG  sys¬ 
tem  mainly  consists  of  three  components:  TEG  module,  exhaust 
channel  and  heat  sink.  Thus,  the  optimization  of  an  ETEG  system 
mainly  has  two  approaches:  one  aims  at  the  development  of 
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Nomenclature 

A 

cross  section  area  of  inlet  and  outlet  (m2) 

Subscripts/superscripts 

B 

fuel  consumption  rate  (kg  h  ') 

1 

normal  vector  1 

E 

electric  field  intensity  vector  (V  m  ') 

2 

normal  vector  2 

h 

heat  transfer  coefficient  (W  m  2  I<  1 ) 

c 

cold  side 

I 

current  (A) 

c,s 

ceramic,  stainless  steel 

J 

current  density  (A  m  2) 

copper 

copper  unit 

k 

turbulence  kinetic  energy  (J  kg-1) 

eff 

effective 

L 

number  of  exhaust  channel 

ext 

exhaust  gas 

N 

number  of  bafflers 

E 

turbulence  kinetic  energy  dissipation 

n 

normal  vector  of  wall 

h 

hot  side 

P 

pressure  (Pa) 

ICE 

internal  combustion  engine 

P 

power  (W) 

i 

direction  of  x  coordinate 

Q 

heat  (W) 

in 

inlet 

S 

energy  source  (W  m  3) 

J 

direction  of  y  coordinate 

T 

temperature  (K) 

k 

turbulence  kinetic  energy 

u 

x-velocity  (ms4) 

l 

loss 

V 

potential  (V) 

n 

n  unit 

V 

source  term  of  Seebeck  potential  (V  m  2) 

n 

net 

V 

velocity  (ms4) 

0 

out 

Ohmic  voltage 
output  voltage 

Greeks  symbols 

ohm 

Ohmic  potential 

a 

Seebeck  coefficient  (V  K  ') 

P 

p  unit 

fi 

angle  between  bafflers  and  bottom  wall  of  exhaust 

red 

reduction  of  fuel 

channel  (°) 

s 

Seebeck  voltage 

g 

turbulent  dissipation  rate  (m  2  s  3) 

sbk 

Seebeck  potential 

n 

efficiency  (%) 

T 

temperature 

2 

thermal  conductivity  (W  m-1  K  ') 

TEG 

thermoelectric  generator 

fl 

effective  turbulent  viscosity  (Pa  s) 

wb 

bottom  wall  of  TEG  modules 

P 

density  (kgm  3) 
electrical  conductivity  (S  m  ’) 

inner  wall  of  exhaust  channel 

high-efficiency  thermoelectric  materials  and  the  design  of  TEG 
structures  (e.g.  segmented  TEG  modules)  [7-14];  and  the  other  is 
mainly  to  enhance  the  heat  transfer  from  heat  source  to  TEG  and 
from  TEG  to  heat  sink  [15-17],  Numerous  experiments  have  been 
conducted  to  investigate  the  key  factors  that  influence  the 
performance  of  ETEG  system  following  these  two  approaches.  For 
example,  new  exhaust  channel  designs  were  proposed.  Lesage 
et  al.  [18]  placed  spiral  and  panel  inserts  in  exhaust  channels  to 
experimentally  investigate  their  effect  on  the  performance  of 
TEG.  In  the  experiment  of  Rezania  [19],  a  heat  exchanger  with  20 


micro-channels  in  the  heat  sink  was  tested  for  the  performance 
improvement  of  TEG. 

To  gain  deep  insights  of  the  transport  phenomena  in  TEG 
systems,  theoretical  and  numerical  models  have  been  developed 
by  researchers.  The  one-dimensional  (1-D)  numerical  model  of 
Chen  et  al.  [20]  considered  the  variable  material  properties,  and 
this  model  was  further  extended  to  a  three-dimensional  (3-D) 
numerical  model  [21],  Besides,  two  approaches  for  developing 
3-D  TEG  models  were  proposed  in  [22],  Except  the  3-D  models 
for  single  TEG  units  [21,22],  Rowe  and  Gao  [23]  developed  a  1-D 


Table  1 

Four  typical  operating  conditions  of  the  diesel  engine  considered  in  this  study  [28], 


Parameters  Case  1 


Power  output  (kW)  258.3 

Fuel  consumption  rate  (kg  h~')  53.22 

Temperature  of  exhaust  gas  (K)  808.15 

Mass  flow  rate  of  exhaust  gas  (kg  s  ’)  0.2981 

Output  temperature  of  jacket  water  (K)  366.15 

Mass  flow  rate  of  jacket  water  (kg  s  ')  2.1 


Case  2 (base) 


235.8 

47.79 

792.15 

0.2752 

365.65 


Case  3 
176.2 
35.43 
747.15 
0.2235 
364.65 


Case  4 
117.7 
23.91 
693.15 
0.1697 
363.75 


of  air  [29], 


Temperature,  T (K)  Density,  p  (kg m  3)  Specific  heat,  cp  (kj  kg  1  K  ’)  Viscosity,  p  (Pa  s)  Thermal  conductivity,  k  (W  m  1  K  1) 

600  0.580  1.051  3.06  x  10  5  0.0466 

800  0.441  1.099  3.65  x  10  5  0.0557 

1000  0.348  1.141  4.24  x  IQ-5  0.0681 
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Exhaust-based  thermoelectric  generator  waste  het 
recovery  system  (cooling  system  not  included) 


Fig.  1.  Schematic  of  the  exhaust-based  thermoelectric  generator  (ETEG)  system  for  a  diesel  engine. 


Design  parameters  of  a  single  TEG  unit. 


Parameters 

p/n  semiconductors 

Ceramic 

Stainless  steel  References 

Seebeck  coefficient  (V  IT1) 

Constant 

±2.2  x  10  4 

_ 

_ 

[15] 

Variable 

±(2224  ±  930.6  *  T  -  0.9905  *  T2)  x  10  5  (300  K  -  500  K,  T  in  K) 

[32] 

Thermal  conductivity  (Wm  1  K  ’) 

1  Constant 

1.5 

350 

14.99  [15,31,32] 

Variable 

(62,605  -277.7  *  T  ±  0.41 31  *  T2)  x  1 0^  (300  K  -  500  K,  T  in  K) 

[32] 

Electrical  conductivity  (S  m~') 

Constant 

1.0  x  105 

5.9  x  107 

[15] 

Variable 

(5112  +  163.4  *  T± 0.6279  *  T2)-’  x  1010  (300  K  -  500  K,  T  in  K) 

[32] 

2x2 

5.2x2 

5.2x2 

Thickness  (mm) 

1.25 

0.5 

0.635 

1 

Exhaust  inlet  (constant  velocity  ) 


Outlet  (surrounding  pressure) 


Coupled  wall  (ceramic  with  exhaust  channel) 


model  to  assess  the  performance  of  TEG  modules  (a  TEG  module 
includes  many  TEG  units).  Similarly,  other  1-D  models  for  TEG 
modules  were  also  developed  [24-26], 


To  investigate  more  detailed  transport  phenomena,  such  as  the 
3-D  fluid  flow  in  heat  source  or  heat  sink  and  the  temperature  dis¬ 
tribution  in  TEG  units,  3-D  numerical  models  are  needed,  and  yet, 
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O  Tj,=3 1  OK,  TC=300K,  Experiment 

—  Th=3 1  OK,  TC=300K,  Model 

□  Tf,=320K,  TC=300K,  Experiment 

-  -  Th=320K,  TC=300K,  Model 
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Fig.  3.  Comparison  of  the  experimental  data  in  [36]  and  simulation  results  [22], 


Simulation  cases. 


Parameters  Specifications 

Cross  section  of  inlet  and  outlet,  A  (mm2)  30  x  40,  60  x  40,  120  x  40 

Number  of  exhaust  channel,  L  1,2,3 

Number  of  bafflers,  N  0, 1,  2, 3, 4,  5,  6 

Angle  of  baffler,  p  (°)  10°,  20°,  30°,  40° 

Operating  conditions  of  ICE  Cases  1-4  of  Table  1 


□  LI 

30  mm  x  40  mm  60  mm  *  40  mm 


Fig.  4.  Different  inlet/outlet  cross  sections  of  exhaust  channel. 


few  works  developed  3-D  numerical  models  for  ETEG  systems. 
Zhou  et  al.  [15]  incorporated  water-fed  heat  exchangers  with 
parallel  and  counter  flow  designs  in  their  TEG  model,  and  dis¬ 
cussed  the  effect  of  heat  exchanger  design  on  performance.  Weng 
et  al.  [27]  developed  a  3-D  exhaust  channel  model  with  a  1-D  TEG 
model  to  predict  the  output  power  of  ETEG  system.  To  the  best  of 
the  authors’  knowledge,  full  3-D  numerical  model  for  ETEG  system 
with  exhaust  channel  and  TEG  modules  were  rarely  developed. 

In  this  study,  a  3-D  numerical  model  for  ETEG  system  of  ICE  is 
developed.  In  the  numerical  simulations,  the  effect  of  the  exhaust 
channel  size  on  the  overall  performance  of  ETEG  system  is  investi¬ 
gated.  Bafflers  are  also  installed  in  the  exhaust  channel  at  different 
locations  and  with  different  orientations,  and  its  influence  on  the 
output  power  of  ETEG  and  pressure  drop  along  channel  under 
different  operating  conditions  is  discussed  in  details. 


Current  (A) 


Current  (A) 

Fig.  5.  Change  of  output  voltage  (a),  power  (b)  and  net  power  (c)  with  current  for 
the  ETEG  systems  with  different  exhaust  inlet/outlet  cross  sections. 


Table  5 

Pressure  drop  and  power  loss  through  the  exhaust  channels  with  three  different  inlet/ 


Cross  section 

drop  Ap  (Pa) 

!mVse 

locity  Inlet  are 

(mm2) 

a  A  Power  loss 

Pi(  W) 

30x40 

1643.0 

97.8 

1064 

171.0 

60  x  40 

296.4 

47.2 

2204 

30.8 

120  x  40 

59.1 

23.2 

4484 

6.1 

2.  Model  development 

The  ICE  considered  in  this  study  is  an  in-line  six-cylinder 
turbocharged  diesel  engine,  and  the  operating  parameters  of  this 


engine  are  given  in  Table  1  [28],  Note  that  the  power  output 
defined  in  Table  1  is  the  actual  power  measured  by  dynamometer. 
The  output  temperature  of  jacket  water  is  measured  at  the  engine 
outlet,  which  is  considered  to  be  the  inlet  condition  of  the  cold  side 
of  the  ETEG  system.  Based  on  the  mass  flow  rate  of  jacket  water 
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Current  (A) 
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Fig.  6.  Change  of  output  voltage  per  exhaust  channel  (a),  power  per  exhaust  channel  (b),  net  power  per  exhaust  channel  (c),  voltage  of  all  exhaust  channel(s)  (d),  power  of  all 
exhaust  channel(s)  (e),  and  net  power  of  all  exhaust  channel(s)  (f)  with  current  for  the  ETEG  systems  with  different  numbers  of  exhaust  channels  per  cylinder  (1  =  1,2,  3). 


Table  6 

Pressure  drop  and  power  loss  through  exhaust  channels  with  different  numbers  of 
exhaust  channels  for  one  cylinder  of  the  six-cylinder  engine. 

Number  Inlet  Pressure  drop  Power  loss  Power  loss  in  all 

of  velocity  per  channel  A p  per  channel  channels  I  •  Pi 

channels  ^(ms-1)  (Pa)  P,(W)  (W) 


1  47.2  296.4  30.8  30.8 

2  22.9  88.7  4.5  8.9 

3  15.3  45.0  1.5  4.6 


and  the  heat  transfer  rate  from  TEG  to  jacket  water,  it  can  be  esti¬ 
mated  that  the  temperature  of  jacket  water  may  increase  by  about 
3-5  K  by  flowing  through  the  ETEG  system  considered  in  this 


study.  After  that,  since  the  jacket  water  still  needs  to  be  cooled 
down  through  radiators  before  entering  the  engine  again,  such 
slight  temperature  increment  of  jacket  water  is  assumed  to  have 
negligible  effect  on  the  engine.  Properties  of  air  are  used  for  the 
engine  exhaust  for  consistent  comparison  of  the  different 
simulation  cases,  as  given  Table  2  [29],  because  the  exact  exhaust 
components  under  different  operating  conditions  are  unknown, 
and  it  was  recommended  that  the  error  associated  with  using  air 
properties  for  exhaust  is  usually  less  than  2%  [30], 

A  schematic  of  the  ETEG  system  (with  empty  exhaust  channel) 
is  shown  in  Fig.  1.  Since  high  flow  rates  in  exhaust  channel  may 
cause  large  pressure  drops,  the  exhaust  channel  shown  in  this 
figure  only  provides  the  flow  path  of  the  exhaust  from  one 
cylinder,  corresponding  to  one-sixth  of  the  exhaust  flow  of  the 
six-cylinder  engine;  and  to  further  decrease  the  pressure  drop, 
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Current  (A) 


Index  of  TEG  modules 


Fig.  8.  Change  of  output  voltage  (a),  power  (b)  and  net  power  (c)  with  current,  and  Seebeck  voltages  of  different  TEG  modules  at  0.1  A  (d)  for  the  ETEG  systems  with  different 
numbers  of  bafflers. 


Table  7 

Ranges  of  X-coordinate  for  the  five  TEG  modules  shown  in  Fig.  2.  Table  8 

-  Pressure  drop  and  power  loss  through  exhaust  channels  with  different  numbers  of 

Index  of  TEG  1  2  3  4  5  bafflers. 


Range  of  X-  -0.3  to  76.7  to  153.7  to  230.7  to  307.7  to 

coordinate  51.7  128.7  205.7  282.7  359.7 

(mm) 


the  cases  with  two  and  three  exhaust  channels  for  one  cylinder  are 
also  investigated.  Details  of  the  flow  rate  and  channel  size  effect 
are  discussed  in  Section  3.  Because  of  the  same  operating 


Number  of  bafflers 


N=1 
N=2 
JV  =  3 


Pressure  drop  Ap  (Pa)  Power  loss  Pi  (W) 


45.0 

58.1 

70.3 

82.4 


113.0 

126.1 

113.4 
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0  50  100  150  200  250  300  350  400 

X  position  (mm) 


(b) 

Solid  line:  temperature  of  wall  of  flow  channel 


X  position  (mm) 


Fig.  9.  Heat  transfer  coefficient  (a)  and  temperatures  of  gas  (middle  of  exhaust  channel  and  1  mm  from  channel  wall)  and  wall  (middle  of  exhaust  channel)  (b)  along  the  flow 
direction  (X-direction)  in  exhaust  channel  with  different  numbers  of  bafflers. 


conditions  of  the  six  cylinders,  only  one  exhaust  channel  with  TEG 
modules  is  considered  in  this  study. 

The  shell  of  the  exhaust  channel  is  stainless  steel  310  with  a 
thickness  of  1  mm.  20  TEG  modules  are  placed  on  the  surface  of 
exhaust  channel  to  be  heated,  and  the  cooling  water  of  the  engine 
is  used  as  the  heat  sink  (not  shown  in  Fig.  1).  A  TEG  module 
includes  160  TEG  units.  These  units  are  electrically  connected  in 
series,  therefore  they  have  the  same  electrical  current.  On  the  other 
hand,  these  units  share  the  same  heat  source  and  heat  sink,  and 
heat  transfers  through  them  simultaneously,  meaning  that  they 
are  thermally  connected  in  parallel.  The  dimensions  of  an  assem¬ 
bled  TEG  module  are  41  mm  x  26  mm  x  3.52  mm.  A  single  TEG 
unit  includes  a  p-type  semiconductor  and  an  n-type  semiconduc¬ 
tor,  and  they  are  connected  with  a  conductor  (e.g.  copper).  The 
insulation  of  the  TEG  units  is  achieved  by  placing  substrates  (e.g. 
ceramic)  on  the  top  and  bottom  sides,  as  shown  in  Fig.  1.  The 
design  parameters  of  a  single  TEG  unit  are  given  in  Table  3 
[15,31,32], 

In  this  study,  the  thermoelectric  material  considered  is  bismuth 
telluride  (Bi2Te3),  as  shown  in  Table  3,  both  the  constant  and  var¬ 
iable  properties  of  the  thermoelectric  (p/n  semiconductors)  mate¬ 
rials  are  considered.  The  numeral  simulations  consume  significant 
computational  time,  and  it  saves  about  80%  of  the  computational 
time  if  constant  properties  are  used  for  the  thermoelectric  materi¬ 
als.  Besides,  the  simulations  in  Sections  3.1 -3.4  only  investigate 


the  exhaust  channel  design  effect  with  the  same  heat  source  and 
heat  sink  conditions,  therefore,  the  constant  material  properties 
are  considered  in  these  sections  under  the  same  engine  operating 
condition  (Case  2  in  Table  1).  For  the  simulations  with  variable 
engine  operating  conditions  (Cases  1-4  in  Table  1)  in  Section  3.5, 
since  the  heat  source  and  heat  sink  change,  the  variable  material 
properties  are  used. 

2.1.  Conservation  equations  for  exhaust  channel 


The  conservation  equations  of  mass,  momentum  and  energy  for 
the  steady-state  flow  in  exhaust  channel  are 

V-t>  =  0  (1) 

V-(vv)  =  -ivp  +  V-(vVv)  (2) 

V-(i>Tex[)  =  V-(Aex[VTex[)  (3) 

In  this  study,  the  Reynolds  number  ranges  from  6000  to  38,000 
in  the  exhaust  channel,  and  the  renormalization  group  (RNG)  k  -  e 
turbulence  model  is  used: 

itSpkUi) = i  (a^§)  +Gk+G^-pE  (4) 
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Fig.  10.  Contours  of  temperature  in  the  middle  planes  of  the  exhaust  channels  (a),  and  near  the  channel  walls  (only  half  wall  is  shown)  with  TEG  modules  (b)  for  different 
numbers  of  bafflers. 


| (p£u,)  -  A  («*, g)  +  c,4 ic,  +  Cs.c.)  -  c^f  - 

(5) 

where  Gk  is  the  turbulence  kinetic  energy  generation  due  to  mean 
velocity  gradient.  Gb  is  the  turbulence  kinetic  energy  generation 
due  to  buoyancy.  ak  and  as  are  the  inverse  effective  Prandtl  numbers 
for  k  and  e.  More  details  of  the  formulation  can  be  found  in  [33], 

2.2.  Conservation  equations  for  thermoelectric  generator 


The  balance  of  electrical  motive  force  due  to  Seebeck  effect  in 
p/n  units  is  written  as 

W*  =  -a p.„VT  (9) 

where  Vsbk  (V)  is  the  Seebeck  potential,  and  v.p  n  (V  K  ')  is  the  See¬ 
beck  coefficient  of  p/n  units.  Taking  the  divergence  of  Vsbk,  Eq.  (9) 
yields: 

V  ■  (-VVsWt)  +  Vsbk  =  0  (10) 

where  Vsbk  (V  m  2)  is  the  source  term: 


The  steady-state  energy  conservation  equation  in  TEG  is 

V  •  [VcopperVTp,n, COpper]  +  ST  =  0  (6) 

where  ipn, capper  (W m  1  K  ')  is  the  thermal  conductivity  of  p-type, 
n-type  and  copper  units.  Tp^copper  (K)  is  the  temperature  of  p-type, 
n-type  and  copper  units.  ST  (W  m  3)  represents  the  energy  source 
such  as  Joule  heat  and  Thomson  heat  generated  in  TEG: 

(~Jp  -  VapTp/p;  p-type  unit 

3-fl  -  n-type  unit  (7) 

wife  flpp.-:  copper  units 


Vsbk 


V  •  [-c/.pVTp] ;  p-type  unit 

V  •  [-a„VT„];  n-type  unit 


The  Ohmic  voltage  drop  due  to  the  current  flow  is 


(11) 


Wokm  =  -Pp;n,coppJp,n,copper  (12) 

Based  on  the  continuity  of  current,  i.e.  V  •  Jp,n,copPer  =  0,  Eq.  (12) 
is  written  as 


V-  copper Wohm]  =  0 


(13) 


2.3.  Boundary  conditions  for  exhaust  channel 


where  Jp,n,copper  (A  m  2)  is  the  electrical  current  intensity  vector,  and 
rr p,n,copper  (S  m  ')  is  the  electrical  conductivity. 

The  energy  conservation  equation  in  ceramic  and  stainless  steel 
310  in  steady  state  only  accounts  for  the  heat  conduction: 

V  •  (4,sVTCiS)  =  0  (8) 

where  (W  m  1  K  ')  is  the  thermal  conductivity  of  ceramic  and 
stainless  steel  units. 


At  the  inlet  of  exhaust  channel,  the  temperature  and  velocity  of 
exhaust  are  defined;  and  at  the  outlet,  the  pressure  is  fixed  at 
1  atm.  Due  to  the  symmetrical  geometry  of  the  exhaust  channel, 
a  quarter  of  the  channel  is  considered  as  the  computational  domain 
with  symmetry  boundary  conditions,  as  shown  in  Fig.  2.  At  the 
gas-solid  interface  in  the  exhaust  channel,  a  coupled  thermal 
boundary  condition,  without  temperature  jump  at  the  interface 
and  satisfying  the  conservation  of  heat  flux,  is  defined: 
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z  Heat  transfer  coefficient  ( W  m'1  KT1) 
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Fig.  11.  Contour  of  heat  transfer  coefficient  for  the  exhaust  channel  walls  (only  half  wall  is  shown)  with  TEG  modules  for  different  numbers  of  bafflers. 


Mt„-r..).(^)_  64) 

where  ha  (W  m  2  K  1 )  is  the  heat  transfer  coefficient  between 
exhaust  and  wall  of  exhaust  channel,  nr  is  the  normal  vector  to 
the  interface.  Text  (K)  is  the  temperature  of  exhaust,  Twi  (K)  is  the 
temperature  at  the  interface,  and  Xs  (W  m  1 1<  ')  is  the  thermal  con¬ 
ductivity  of  the  channel  shell.  The  heat  loss  of  the  exhaust  through 
the  other  channel  walls  is  neglected  in  this  study.  At  the  interface 
between  exhaust  channel  shell  and  TEG  module,  it  also  satisfies  that 
the  temperate  is  continuous  and  the  heat  flux  is  conserved. 

2.4.  Boundary  conditions  for  thermoelectric  generator 

A  constant  heat  transfer  coefficient  (hc,  1 000  W  m  2  I<  1 )  is 
defined  on  the  cold  side  of  TEG  modules  representing  the  engine 
coolant  cooling  effect,  which  is  estimated  using  the  Gnielinski 
equation  [34]  and  was  also  used  in  previous  TEG  models  (e.g. 
[35]).  The  heat  flux  to  the  heat  sink  is  calculated  as 

(,5) 

where  hc  (W  m  2  K  ’)  is  the  heat  transfer  coefficient  between  heat 
sink  and  cold  side  of  TEG  modules;  n2  is  the  normal  vector  of  the 
interface  of  ceramic  and  heat  sink,  which  points  at  heat  sink;  Twb 


(K)  is  the  wall  temperature  of  cold  side  of  TEG  modules;  Tc  (K)  is 
the  temperature  of  heat  sink;  and  Xc  (W  m  1  K  1 )  is  the  thermal 
conductivity  of  cold  side  of  TEG  modules,  i.e.  ceramic. 

As  mentioned  previously,  a  TEG  module  includes  a  number  of 
single  TEG  units,  and  the  TEG  units  are  connected  thermally  in 
parallel  and  electrically  in  series.  Therefore,  all  the  TEG  units  have 
the  same  operating  current.  In  each  TEG  unit,  for  the  Ohmic  poten¬ 
tial,  the  operating  current  density  and  a  reference  potential  of  zero 
are  defined  at  the  two  ends  of  the  copper  conductors,  respectively; 
for  the  Seebeck  potential,  an  electric  field  intensity  vector  and  a 
reference  potential  of  zero  are  defined  at  the  two  ends  of  the  p 
or  n  units,  respectively.  The  electric  field  intensity  vector,  £ 
(V  m_1),  is  defined  as 

£=ap,„VTp,n  (16) 

The  details  of  the  single  TEG  unit  model  can  also  be  found  in 
[22],  The  output  voltage  of  a  TEG  module  is  the  summation  of 
the  Seebeck  voltage  and  the  Ohmic  voltage  drop  of  all  TEG  units 
in  the  module: 

V„ut  =  5>  +  £V„  (17) 

where  Vout  (V)  is  the  output  voltage  of  a  TEG  module,  Vs  (V)  is  the 
Seebeck  voltage  of  a  single  TEG  unit,  and  Va  (V)  is  the  Ohmic  voltage 
drop  of  a  single  TEG  unit. 
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2.5.  Numerical  methods 

This  numerical  model  is  implemented  in  the  computational 
fluid  dynamics  software  FLUENT  by  using  its  user  defined  functions 
to  customize  the  formulation.  The  individual  variables  of  the  con¬ 
servation  equations  are  solved  by  using  pressure-based  segregated 
solver.  The  second  order  upwind  method  is  used  for  spatial  discret¬ 
ization.  An  algebraic  multigrid  (AMG)  method  with  a  Gauss-Seidel 
type  smoother  is  used  to  accelerate  the  convergence.  The  compu¬ 
tational  domain  and  mesh  with  some  of  the  boundary  conditions 
are  illustrated  in  Fig.  2.  The  grid  independency  study  carried  out 
indicates  that  4,118,400  cells  for  five  TEG  modules  and  2,059,680 
cells  for  one  exhaust  channel  are  sufficient. 

3.  Results  and  discussion 

The  TEG  modeling  results  have  been  compared  with  the  exper¬ 
imental  data  in  [36]  previously  [22],  with  reasonable  agreement,  as 
shown  in  Fig.  3.  It  should  be  noticed  that  since  the  current  changes 
almost  linearly  with  voltage,  the  slope  of  the  curve  is  mainly  deter¬ 
mined  by  the  electrical  resistance,  and  the  magnitude  is  mainly 
determined  by  the  temperature  difference  (mainly  depends  on 
the  thermal  resistance)  and  Seebeck  coefficient.  Therefore,  with 
fixed  material  properties,  by  properly  considering  the  electrical 
and  thermal  contact  resistances  at  the  interfaces  between  different 
materials,  it  is  possible  to  have  good  agreement  between  modeling 
and  experimental  results.  In  this  study,  the  exhaust  channel  is  con¬ 
sidered  with  a  large  number  of  TEGs  in  the  simulations  represent¬ 
ing  practical  3-D  ETEG  system.  The  simulation  cases  are 
summarized  in  Table  4.  Three  cross-section  areas  of  inlet  and  out¬ 
let,  A  (m2),  are  compared  first.  Then,  the  effect  of  number  of 
exhaust  channels  corresponding  to  one  cylinder  of  the  six-cylinder 


engine,  L,  is  tested.  After  that,  bafflers  are  mounted  in  the  exhaust 
channel  to  enhance  the  performance  of  ETEG.  The  influence  of 
number  of  bafflers,  N,  is  investigated,  and  an  improved  exhaust 
channel  design  is  proposed.  The  cases  with  different  angles 
between  baffler  and  inner  wall  of  exhaust  channel,  /?  (°),  are  also 
simulated  for  design  optimization.  Finally,  to  assess  the  perfor¬ 
mance  of  ETEG  under  the  different  operating  conditions  of  the 
ICE,  the  four  operating  conditions  given  in  Table  1  are  simulated 
and  compared. 

3.1.  Effect  of  exhaust  channel  inlet/outlet  area 

Three  inlet/outlet  cross-section  areas  of  exhaust  channel  are 
compared  in  this  section,  as  shown  in  Fig.  4  and  Table  4.  The  only 
difference  among  the  three  cross  sections  is  the  width,  and  for 


Fig.  13.  Schematic  of  different  baffler  angles  (/I). 
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consistent  comparison,  all  the  other  design  and  operating 
parameters  are  kept  the  same.  For  this  section,  one  exhaust  chan¬ 
nel  is  considered  for  one  cylinder  of  the  six-cylinder  engine. 

The  voltage,  power  and  net  power  outputs  at  different  currents 
of  the  TEGs  are  shown  in  Fig.  5a-c.  It  should  be  noticed  that  the 
current-voltage  and  current-power  curves  in  this  study  are  all 
obtained  from  numerical  simulations  in  steady-state  operating 
conditions  and  with  the  same  numerical  methods,  therefore  the 
numerical  results  are  all  repeatable.  The  statistical  error  is 
neglected  in  the  numerical  simulations,  which  is  also  common  in 
previous  numerical  models  (e.g.  [27,28]).  This  is  different  from 
experimental  measurements,  in  which  the  measured  values  in  dif¬ 
ferent  experiments  may  change,  and  the  statistical  error  needs  to 
be  considered.  For  all  the  three  designs,  the  output  voltage 
decreases  almost  linearly  with  the  increment  of  current,  and  the 
power  and  net  power  follow  the  increasing-decreasing  trend.  Note 
that  the  net  power  is  the  TEG  output  power  minus  the  power  loss 
of  exhaust  through  the  channel.  The  power  loss  Pi  (W)  caused  by 
the  pressure  drop  A P  (Pa)  from  inlet  to  outlet  is  calculated  as 
P[  =  Ap  A  -  vin  (18) 

Pn=Po~Pl  (19) 

where  A  (m2)  is  the  cross-section  area  of  inlet/outlet,  vin  (m  s  ’)  is 
the  velocity  of  exhaust  at  inlet,  P„  (W)  is  the  net  output  power, 
and  P0  (W)  is  the  TEG  output  power. 

The  pressure  or  power  loss  in  exhaust  channel  may  affect  the 
engine  operation,  and  the  recovery  of  such  power  loss  requires 
additional  pumping  power.  Therefore,  the  net  power  output  could 
better  reflect  the  overall  ETEG  performance  than  the  TEG  power 
output. 


With  the  smallest  cross  section  of  30  mm  x  40  mm,  the  voltage 
and  power  are  the  highest  (Fig.  5a  and  b).  Taking  3  A  as  an  exam¬ 
ple,  the  output  voltages  the  TEGs  with  cross  sections  of 
30  mm  x  40  mm,  60  mm  x  40  mm  and  120  mm  x  40  mm  are 
56.1  V,  50.6  V  and  31.5  V,  respectively.  The  reason  is  that  with 
the  same  exhaust  mass  flow  rate,  a  smaller  cross  section  increases 
the  flow  velocity,  leading  to  enhanced  heat  transfer  to  the  TEGs. 

Flowever,  in  considering  the  power  loss  and  pressure  drop 
through  exhaust  channel  (Table  5),  the  exhaust  channel  with  the 
smallest  cross  section  of  30  mm  x  40  mm  has  the  largest  pressure 
drop,  about  six  times  that  with  the  cross  section  of 
60  mm  x  40  mm,  and  28  times  that  of  the  cross  section  of 
120  mm  x  40  mm,  due  to  the  high  flow  velocity  in  the  small  cross 
sections.  As  a  result,  the  net  power  of  the  ETEG  with  the  smallest 
cross  section  of  30  mm  x  40  mm  is  the  lowest,  and  the  TEG  output 
power  cannot  even  overcome  the  power  loss  through  the  exhaust 
channel,  as  shown  in  Fig.  5c.  The  moderate  cross  section  of 
60  mm  x  40  mm  produces  the  highest  net  power  output  among 
the  three  designs  for  most  of  the  operating  currents,  and  only  when 


Table  9 

Pressure  drop  and  ] 


Baffler  angle 
Empty 
0  =  10° 

0  =  20° 

0  =  30° 

0  =  40° 


loss  through  exhaust  channels  with  different  baffler  angles. 


Pressure 

45.0 

67.7 

113.4 

251.8 

675.2 


drop  Ap  (Pa) 


Power  loss  P,  (W) 

1.5 
2.3 
3.8 

8.5 
22.8 


Fig.  14.  Change  of  output  voltage  (a),  power  (b)  and  net  power  (c)  wit 
baffler  angles. 


t,  and  Seebeck  voltages  of  different  TEG  modules  at  0.1  A  (d)  1 
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X  position  (mm) 


X  position  (mm) 


(C) 


X  position  (mm) 


Fig.  15.  Heat  transfer  coefficient  (a),  gas  temperature  (middle  of  exhaust  channel 
and  1  mm  from  channel  wall)  (b),  and  wall  temperature  (middle  of  exhaust 
channel)  (c)  along  the  flow  direction  (X-direction)  in  exhaust  channel  with  different 
baffler  angles. 


the  current  is  low  (less  than  1.3  A),  the  largest  cross  section  of 
120  mm  x  40  mm  produces  the  highest  net  output  power. 

Therefore,  the  exhaust  channel  with  the  cross  section  of 
60  mm  x  40  mm  might  be  considered  to  balance  the  flow  resis¬ 
tance  of  the  TEG  power  output,  and  this  cross  section  is  used  for 
the  rest  of  the  simulations  in  this  study. 

3.2.  Effect  of  number  of  exhaust  channels 

As  mentioned  previously,  the  ICE  considered  in  this  study  is  a 
six-cylinder  diesel  engine.  Since  high  flow  rates  in  exhaust  channel 
may  cause  large  pressure  drops,  in  this  study,  one,  two  and  three 
exhaust  channels  (1  =  1,  2,  3)  are  considered  to  only  provide  the 


flow  path  for  the  exhaust  from  one  cylinder,  corresponding  to 
one-sixth  of  the  exhaust  flow  of  the  six-cylinder  engine. 

Based  on  the  base  operating  condition  (Case  2  in  Table  1 ),  the 
cases  for  one,  two  and  three  exhaust  channels  (L  =  1,  2,  3)  to  pro¬ 
vide  the  flow  path  of  the  exhaust  from  one  cylinder  are  studied. 
In  fact,  the  only  difference  of  the  simulation  cases  with  different 
numbers  of  exhaust  channels  (I)  is  the  velocity  at  the  inlet  of 
exhaust  channel.  Fig.  6  shows  the  change  of  voltage,  power  and 
net  power  with  current  for  each  exhaust  channel  and  for  all  the 
exhaust  channels.  If  the  average  voltage/power  of  the  exhaust 
channels  for  the  cases  with  two  and  three  channels  are  used  for 
comparison  (Fig.  6a-c),  the  single-exhaust  channel  case  shows 
both  the  highest  power  and  net  power,  and  they  decrease  with 
more  exhaust  channels.  For  example,  at  3  A,  the  average  voltages 
are  50.6  V,  28.7  V  and  16.3  V  for  one,  two  and  three  channels, 
respectively.  The  average  voltage  of  the  single-channel  case  is 
about  3  times  of  the  three-channel  case.  With  more  exhaust  chan¬ 
nels,  the  peak  power  is  produced  at  lower  current,  suggesting  that 
increasing  the  number  of  channels  reduces  the  performance  of 
each  TEG.  However,  in  considering  the  power  loss  in  exhaust  chan¬ 
nel  in  Table  6,  it  is  much  lower  for  the  three-channel  case  than  the 
single-channel. 

As  a  result,  Fig.  6d-f  shows  that  the  three-channel  case  pro¬ 
duces  both  the  highest  TEG  power  and  net  power,  especially  for 
the  net  power.  It  can  also  be  noticed  that  the  improvement  in  total 
net  power  by  increasing  the  number  of  exhaust  channels  becomes 
less  significant  with  more  channels.  In  addition,  with  more  exhaust 
channels,  the  current  corresponding  to  peak  power  decreases. 
Therefore,  properly  increasing  the  number  of  exhaust  channels 
may  improve  the  ETEG  performance,  however,  since  more  space 
and  TEG  modules  are  needed,  the  system  size  and  cost  need  to 
be  considered  as  well.  It  is  also  suggested  that  corresponding  to 
different  engine  operating  conditions  (exhaust  flow  rates),  making 
the  number  of  channels  adjustable  may  further  improve  the  ETEG 
performance. 

The  results  in  Sections  3.1  and  3.2  show  that  the  size  and  shape 
of  exhaust  channel  may  affect  ETEG  performance  significantly,  and 
proper  channel  design  is  of  paramount  importance,  because  the 
heat  transfer  rate  from  exhaust  to  TEG  and  the  pressure  drop 
through  exhaust  channel  change  significantly  with  channel  design. 

3.3.  Effect  of  number  of  bafflers 

In  this  section,  the  influence  of  number  of  bafflers  mounted  in 
exhaust  channel  is  investigated.  The  bafflers  are  mounted  to 
enhance  the  heat  transfer  to  the  TEG  modules.  The  number  of 
baffler,  defined  as  N,  ranges  from  0  to  6  per  channel.  The  dimension 
of  each  baffler  is  1 58  mm  x  20  mm  x  1  mm,  along  the  X-,  Y-  and 
Z-directions,  respectively.  The  angle  between  the  channel  wall 
and  the  baffler,  ft,  is  20°.  All  the  bafflers  are  evenly  spaced  in  the 
exhaust  channel.  A  detailed  schematic  of  the  exhaust  channel  with 
different  numbers  of  bafflers  is  shown  in  Fig.  7.  The  empty  exhaust 
channel  without  baffler  is  considered  as  the  base  case  for  compar¬ 
ison.  For  consistent  comparison,  except  the  bafflers,  all  the  other 
design  and  operating  parameters  are  kept  the  same  for  the  simula¬ 
tions  in  this  section. 

It  can  be  observed  in  Fig.  8a-c  that  the  voltage  decreases  almost 
linearly  as  the  current  increases  for  all  the  cases,  and  the  power 
and  net  power  all  follow  the  increasing-decreasing  trend.  The  peak 
output  power  is  obtained  at  about  2  A  for  all  the  cases.  With  one 
baffler  mounted  at  the  inlet,  the  net  power  is  the  lowest,  and  only 
when  the  number  of  bafflers  becomes  more  than  3,  the  net  power 
is  higher  than  the  empty  channel  case.  As  shown  in  Fig.  8d,  for  the 
Seebeck  voltages  of  the  five  TEG  modules,  with  one  baffler  at  the 
inlet,  the  Seebeck  voltage  of  TEG  1  is  increased,  however,  it  is  even 
lower  than  the  empty  channel  case  for  TEGs  2-5.  With  more 
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bafflers  mounted,  the  Seebeck  voltages  of  the  TEG  modules  down¬ 
stream  are  increased. 

The  ranges  of  X-coordinate  of  the  five  TEG  modules  (Fig.  2)  are 
given  in  Table  7,  showing  that  TEGs  1-5  are  mounted  along  the 
flow  direction.  Although  more  bafflers  improve  the  performance, 
the  pressure  drop  along  the  exhaust  channel  increases  as  well,  as 
shown  in  Table  8  with  an  almost  linear  relationship.  Since  the 
results  in  Fig.  8  suggest  that  bafflers  need  to  be  placed  in  the  whole 
channel  to  ensure  the  performance  of  all  the  TEG  modules,  using 
five  bafflers  rather  than  six  corresponding  to  the  five  TEG  modules 
rather  than  six,  and  placing  the  five  bafflers  just  before  the  five  TEG 
modules  (the  “N  =  5  im”  case)  shows  similar  net  power  (Fig.  8)  and 
lower  pressure  drop  (Table  8)  than  the  six-baffler  case. 

To  understand  the  effect  of  number  of  bafflers  shown  in  Fig.  8, 
the  heat  transfer  coefficient  (defined  in  Eq.  (14))  and  the  tempera¬ 
tures  of  gas  near  channel  wall  and  wall  along  the  flow  direction 
(X-direction)  in  exhaust  channel  with  different  numbers  of  bafflers 
are  shown  in  Fig.  9.  The  exhaust  gas  temperature  shown  in  this 
figure  is  in  the  middle  of  the  channel,  along  the  line  1  mm  away  from 
the  channel  wall  with  the  TEG  modules,  and  the  wall  temperature  is 
along  the  middle  line  of  the  wall  with  the  TEG  modules.  It  can  be 
noticed  that  the  heat  transfer  coefficient  (Fig.  9a)  is  greatly  enhanced 
at  the  locations  with  bafflers,  and  keeps  almost  constant  at  the 
locations  without  bafflers.  Even  only  with  one  baffler  at  the  inlet, 
the  heat  transfer  coefficient  is  enhanced  for  the  whole  channel  com¬ 
pared  with  the  empty  channel  case.  However,  the  near  wall  gas 
temperature  (Fig.  9b)  decreases  more  significantly  than  the  empty 
channel  case  for  the  whole  channel.  As  a  result,  compared  with  the 
empty  channel  case,  the  near  wall  exhaust  gas  temperature  is  only 
higher  at  the  first  TEG  module,  but  lower  downstream.  Correspond¬ 
ingly,  the  wall  temperature  changes  with  the  near  wall  exhaust  gas 
temperature,  influencing  the  performance  of  the  TEG  modules. 

The  reason  is  that  with  only  one  baffler  at  the  inlet,  although  the 
flow  velocity  is  enhanced  near  the  wall,  only  part  of  the  hot  exhaust 
gas  is  guided  to  the  near  wall  region.  As  a  result,  after  the  heat  is 
absorbed  by  the  first  TEG  module,  the  near  wall  exhaust  gas  temper¬ 
ature  becomes  much  lower  than  the  exhaust  gas  temperature  in  the 
middle  of  the  channel.  As  shown  in  Fig.  1  Oa,  the  empty  channel  case 


shows  more  evenly  distributed  temperature  than  the  one  baffler 
case,  and  with  more  bafflers  downstream,  more  hot  exhaust  gases 
are  forced  to  flow  to  the  near  wall  region,  improving  both  the  heat 
transfer  coefficient  and  near  wall  exhaust  gas  temperature.  Corre¬ 
spondingly,  the  near  wall  temperature  shown  in  Fig.  10b  shows  the 
similar  trend.  Moreover,  since  the  bafflers  installed  fully  covers  the 
exhaust  channel  along  the  Z-direction,  the  heat  transfer  coefficient 
only  changes  along  the  flow  direction  (X-direction),  and  the  change 
along  the  Z-direction  is  insignificant  (Fig.  1 1 ).  In  addition,  the  velocity 
magnitude  contours  shown  in  Fig.  1 2  further  confirm  the  baffler  effect 
on  the  gas  flow  in  exhaust  channel. 

The  results  in  this  section  suggest  that  to  enhance  the  heat 
transfer  from  exhaust  channel  to  TEG  modules,  bafflers  need  to 
be  placed  at  all  the  locations  with  TEG  modules.  Although  only 
placing  bafflers  at  the  channel  inlet  could  increase  the  heat  transfer 
coefficient  for  the  whole  channel,  the  near  wall  temperature  down¬ 
stream  might  decrease  significantly,  leading  to  performance  degra¬ 
dation  of  the  TEG  modules  downstream.  If  bafflers  can  be  placed 
properly  according  to  the  locations  of  the  TEG  modules,  further 
increasing  the  number  of  bafflers  may  not  improve  the  TEG  perfor¬ 
mance  significantly,  but  further  increase  the  pressure  drop  along 
channel.  For  example,  for  the  ETEG  system  considered  in  this  study, 
it  is  suggested  to  place  one  baffler  just  before  each  TEG  module. 

3.4.  Effect  of  baffler  angle 

In  this  section,  the  effect  of  the  angle  between  baffler  and  bot¬ 
tom  wall  of  exhaust  channel  is  investigated,  which  is  defined  as 
p,  as  shown  in  Fig.  13.  Based  on  the  improved  exhaust  channel  of 
“N=  5  im”  obtained  in  Section  3.3,  four  kinds  of  exhaust  channels 
with  different  baffler  angles  are  considered  (Fig.  13).  The  empty 
exhaust  channel  without  baffler  is  considered  as  well  for 
comparison.  For  consistent  comparison,  except  the  bafflers,  all 
the  other  design  and  operating  parameters  are  kept  the  same  for 
the  simulations  in  this  section. 

An  obvious  improvement  of  the  output  voltage  and  power  is 
observed  as  the  baffler  angle  increases,  as  shown  in  Fig.  14a  and 
b.  Taking  the  current  of  0.1  A  as  an  example,  the  output  voltages 
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Fig.  17.  Contours  of  temperature  in  the  middle  planes  of  the  exhaust  channels  (a),  and  near  the  channel  walls  (only  half  wall  is  shown)  with  TEG  modules  (b)  for  different 
baffler  angles. 


of  the  empty  channel,  baffler  angles  of  10°,  20°,  30°  and  40°  are 
69.4  V,  69.7  V,  75.2  V,  82.4  V  and  92.8  V,  respectively.  The  percent¬ 
ages  of  improvement  compared  with  the  empty  channel  case  are 
0.4%,  8.4%,  18.7%  and  33.7%,  respectively.  Although  Table  9  shows 
that  the  pressure  drop  is  significantly  increased  with  the  increment 
of  baffler  angle,  especially  from  30°  to  40°,  the  net  power  is  still  the 
highest  for  the  angle  of  40°,  as  shown  in  Fig.  14c.  With  the  smallest 
angle  of  10°  shown  in  Fig.  14d,  the  Seebeck  voltages  at  0.1  A  of  the 
first  three  TEG  modules  are  improved  compared  with  the  empty 
channel  case,  however,  the  last  two  TEG  modules  have  lower  See¬ 
beck  voltages  with  bafflers.  The  reason  is  that  with  enhanced  heat 
transfer  at  the  inlet  of  channel,  more  downstream  hot  gas  needs  to 
be  guided  to  flow  to  the  near  wall  region,  and  the  angle  of  10°  is  not 
sufficient  for  that  purpose.  The  heat  transfer  coefficient  and 
temperature  shown  in  Fig.  15  confirm  that  although  heat  transfer 
coefficient  is  improved  even  with  the  baffler  angle  of  10°,  the  near 
wall  gas  temperature  downstream  becomes  lower  than  the  empty 
channel  case,  suggesting  that  the  baffler  angle  needs  to  be 


sufficiently  high  to  ensure  the  near  wall  gas  temperature.  With  a 
higher  baffler  angle,  near  wall  gas  temperature  variation  along 
the  channel  becomes  more  significant,  with  the  trend  that  the 
temperature  decreases  due  to  the  TEG  module  heat  absorption, 
and  increases  because  the  baffler  guide  more  hot  gas  flowing  to 
the  near  wall  region. 

Fig.  16  shows  that  the  heat  transfer  coefficient  is  almost  evenly 
distributed  for  the  empty  channel  case,  and  with  the  increment  of 
baffler  angle,  the  heat  transfer  coefficient  becomes  more  unevenly 
distributed.  Compared  by  the  empty  channel  case,  Fig.  17  shows 
that  the  baffler  angle  of  10°  results  in  more  evenly  distributed 
temperature  along  the  Y-direction,  and  by  further  increasing  the 
baffler  angle,  more  hot  gas  in  the  middle  is  able  to  flow  to  the  near 
wall  region,  the  temperature  becomes  more  evenly  distributed 
along  the  Y-direction  and  higher  near  the  wall,  leading  to  improved 
performance  of  the  TEG  modules.  The  velocity  magnitude  distribu¬ 
tion  shown  in  Fig.  18  confirms  the  improvement  of  heat  transfer 
coefficient  with  larger  baffler  angles.  With  larger  baffler  angles,  it 
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Fig.  18.  Contours  of  velocity  magnitude  in  the  middle  planes  of  the  exhaust  channels  (a),  and  near  the  channel  walls  (only  half  wall  is  shown)  with  TEG  modules  (b)  for 
different  baffler  angles. 


can  also  be  noticed  that  the  velocity  magnitude  is  significantly 
increased,  which  is  the  reason  of  the  high  pressure  drop  in  Table  9. 

To  further  assess  the  performance  of  the  ETEG  system,  two 
kinds  of  efficiency  are  calculated.  One  is  the  efficiency  of  TEG: 

=  (20) 


where  Qjeg  (W)  is  the  heat  transfer  rate  into  the  hot  side  of  TEG,  and 
Pn  (kW)  is  the  net  power  output  of  ETEG  system.  Another  one  is  the 
efficiency  of  exhaust  waste  heat  recovery: 


where  CW  (W)  is  the  heat  removal  rate  from  engine  carried  by 
exhaust  gas  (only  part  of  this  heat  can  transfer  into  the  hot  side 
of  TEG).  The  reduction  of  engine  fuel  consumption  by  using  ETEG 
system  can  also  be  estimated: 


Bred  =  (22) 

where  Pice  (kW)  is  the  engine  power  output  (Table  1),  and  Bice 
(kg  h  ')  is  the  engine  fuel  consumption  rate  (Table  1). 

Fig.  19  shows  the  reduction  of  fuel  consumption  rate  per  cylin¬ 
der,  waste  heat  recovery  efficiency  of  exhaust  gas,  and  efficiency  of 
TEG  modules  at  different  operating  currents  of  the  TEG  modules 
with  the  baffler  angle  (/i)  of  30°.  Taking  2  A  as  an  optimal  example, 
the  reduction  of  fuel  consumption  rate  per  cylinder,  waste  heat 
recovery  efficiency  of  exhaust  gas  and  efficiency  of  TEG  modules 
are  about  50 g h  ',  1%  and  4%,  respectively.  For  the  six-cylinder 
engine,  the  total  reduction  of  fuel  consumption  rate  is  therefore 
about  0.3  kg  h  ',  meaning  that  the  fuel  consumption  rate  of  the 
engine  is  reduced  by  about  0.6%  by  using  the  ETEG  system  for 
waste  heat  recovery.  By  comparing  the  waste  heat  recovery  effi¬ 
ciency  of  exhaust  gas  (1%)  and  efficiency  of  TEG  modules  (4%),  it 
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Fig.  19.  Reduction  of  fuel  consumption  rate  per  cylinder,  waste  heat  recovery 
efficiency  of  exhaust  gas,  and  efficiency  of  TEG  modules  at  different  operating 
currents  of  the  TEG  modules  with  the  baffler  angle  ( ji )  of  30°. 


can  be  found  that  only  about  25%  of  the  waste  heat  in  exhaust  gas 
is  transferred  into  the  TEG  modules,  suggesting  that  there  is  still  a 
significant  potential  to  improve  the  performance  of  ETEG  system 
by  heat  transfer  enhancement. 

The  results  in  this  section  suggest  that  to  ensure  effective 
utilization  of  all  the  hot  exhaust  gas,  the  baffler  angle  needs  to 
be  sufficiently  high,  especially  for  the  downstream  locations.  How¬ 
ever,  higher  baffler  angles  increase  the  pressure  drop  significantly, 
therefore,  it  is  suggested  that  variable  baffler  angles,  with  the  angle 
increasing  along  the  flow  direction,  might  be  a  middle  course  for 
balancing  the  heat  transfer  and  pressure  drop. 

3.5.  Effect  of  engine  operating  condition 

In  this  section,  the  effect  of  the  different  engine  operating  con¬ 
ditions  shown  in  Table  1  is  investigated.  According  to  the  results  in 
Section  3.4,  the  baffler  angle  of  30°  is  considered  in  this  section, 
because  this  angle  could  better  balance  the  heat  transfer  and 
pressure  drop.  With  the  different  engine  operating  conditions, 
the  variable  material  properties  of  the  p/n  semiconductors  in 
Table  3  are  used  in  the  simulations.  The  four  operating  conditions 
in  Table  1  all  correspond  to  the  engine  speed  of  1500  rpm,  and  the 
main  difference  is  the  power  output.  From  Cases  1  to  4,  the  output 
power  is  decreased  from  258.3  kW  to  117.7  kW.  For  consistent 
comparison,  except  the  engine  operating  condition,  all  the  other 
design  and  operating  parameters  are  kept  the  same  for  the  simula¬ 
tions  in  this  section.  Fig.  20a  and  b  shows  that  the  output  voltage 
and  power  of  the  ETEG  system  decrease  from  Case  1  to  Case  4, 
i.e.  decease  with  increment  of  engine  power.  Taking  the  current 
of  0.5  A  as  an  example,  the  output  powers  for  Cases  1-4  are 
23.9  W,  22.4  W,  18.3  W  and  13.1  W,  respectively.  The  main  reason 
is  the  higher  exhaust  gas  temperature  and  flow  rate  with  higher 
engine  output  power,  as  shown  in  Table  1 .  However,  in  considering 
the  pressure  drop  in  exhaust  channel  (Table  10),  Case  2  produces 
the  highest  net  power.  The  results  suggest  that  a  single  ETEG 
design  may  not  be  suitable  to  the  various  engine  operating 
conditions,  and  making  the  number  of  exhaust  channels  and  baffler 
angle  adjustable  according  to  different  engine  operating  condi¬ 
tions,  as  suggested  in  Sections  3.2  and  3.4,  might  be  able  to  further 
improve  the  performance. 


Table  10 

Pressure  drop  and  power  loss  through  exhaust  channels  with  different  engine 
operating  conditions  (Table  1 ). 

Case  Inlet  velocity  vin  (m  s-,3  Pressure  drop  Ap  (Pa)  Power  loss  Pj  (W) 

1  17.4  327.4  12.6 

2  15.3  251.8  8.5 

3  11.9  165.7  4.4 

4  8.4  90.3  1.7 


4.  Conclusion 

In  this  study,  a  3-D  numerical  model  for  exhaust-based  thermo¬ 
electric  generator  (ETEG)  system  has  been  developed.  The  effects  of 
exhaust  channel  size,  number  of  channel,  number  of  baffler,  baffler 
angle  and  engine  operating  condition  on  the  heat,  mass  and 
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electric  transfer  characteristics  and  performance  of  ETEG  are 
investigated  in  details.  By  considering  the  detailed  geometry  of 
the  thermoelectric  generator  (TEG)  and  exhaust  channel  as  the 
computational  domain,  various  transport  phenomena  are  observed 
and  elucidated.  Suggestions  on  design  optimization  are  given 
based  on  the  numerical  results.  It  is  found  that  a  smaller  inlet  of 
exhaust  channel  enhances  the  heat  transfer  to  TEG  modules,  but 
increases  the  flow  resistance,  and  therefore  a  moderate  channel 
size  is  needed.  Based  on  the  design  and  operating  conditions  in  this 
study,  the  exhaust  channel  with  the  cross  section  of 
60  mm  x  40  mm  might  be  optimal  to  balance  the  flow  resistance 
and  TEG  power  output.  Increasing  the  number  of  exhaust  channels 
may  improve  the  ETEG  performance,  however,  since  more  space 
and  TEG  modules  are  needed,  the  system  size  and  cost  need  to 
be  considered  as  well.  With  multiple  channels  per  cylinder,  it  is 
also  recommended  to  make  the  number  of  working  channels 
adjustable  for  further  performance  improvement.  To  enhance  the 
heat  transfer  from  exhaust  channel  to  TEG  modules,  bafflers  need 
to  be  placed  near  all  the  TEG  modules.  Although  only  placing 
bafflers  at  the  channel  inlet  could  increase  the  heat  transfer 
coefficient  for  the  whole  channel,  the  near  wall  temperature  down¬ 
stream  might  decrease  significantly,  leading  to  performance  degra¬ 
dation  of  the  TEG  modules  downstream.  If  bafflers  can  be  placed 
properly  according  to  the  locations  of  the  TEG  modules,  further 
increasing  the  number  of  bafflers  may  not  improve  the  TEG  perfor¬ 
mance  significantly,  but  increase  the  pressure  drop  along  channel. 
For  example,  for  the  ETEG  system  considered  in  this  study,  it  is 
suggested  to  place  one  baffler  just  before  each  TEG  module.  To 
ensure  effective  utilization  of  the  hot  exhaust  gas,  the  baffler  angle 
needs  to  be  sufficiently  large,  especially  for  the  downstream 
locations.  However,  larger  baffler  angles  increase  the  pressure  drop 
significantly,  therefore,  it  is  suggested  that  variable  baffler  angles, 
with  the  angle  increasing  along  the  flow  direction,  might  be  a 
middle  course  for  balancing  the  heat  transfer  and  pressure  drop. 
For  the  six-cylinder  engine  and  ETEG  system  in  this  study,  the  total 
reduction  of  fuel  consumption  rate  can  reach  about  0.3  kg  h  \ 
meaning  that  the  fuel  consumption  rate  of  the  engine  is  reduced 
by  about  0.6%  by  using  the  ETEG  system  for  waste  heat  recovery. 
By  comparing  the  waste  heat  recovery  efficiency  of  exhaust  gas 
(1%)  and  efficiency  of  TEG  modules  (4%),  it  can  be  found  that  only 
about  25%  of  the  waste  heat  in  exhaust  gas  is  transferred  into  the 
TEG  modules,  suggesting  that  there  is  still  a  significant  potential  to 
improve  the  performance  of  ETEG  system  by  heat  transfer 
enhancement.  A  single  ETEG  design  may  not  be  suitable  to  the 
various  engine  operating  conditions,  and  making  the  number  of 
exhaust  channels  and  baffler  angle  adjustable  according  to  differ¬ 
ent  engine  operating  conditions  might  be  able  to  further  improve 
the  performance. 
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